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Abstract. We use cosmological magneto-hydrodynamic simulations to study the evolution of magnetic fields in galaxy clusters 
in two different cosmological models, a standard-CDM and a A-CDM model. We show that the magnetic field strength profiles 
closely follow the cluster density profiles outside a core region of radius ~ 200h~ 1 kpc. The magnetic field has a correlation 
length of order 50/i _1 kpc and reverses on scales of ~ 100/i -1 kpc along typical lines-of-sight. The power spectrum of the 
magnetic field can well be approximated by a steep power law with an exponent of ~ —2.7. The mean magnetic field in the 
cluster cores grows roughly exponentially with decreasing redshift, B ~ 10 _2,5z fj,G. Merger events have a pronounced effect 
on magnetic field evolution, which is strongly reflected in measurable quantities like the Faraday rotation. The field evolution 
in the two different cosmologies proceeds virtually identically. All our cluster models very well reproduce observed Faraday 
rotation measurements when starting with nG seed fields. 
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1. Introduction 

While Faraday -rotation measurements clearly show that galaxy 
clusters contain magnetic fields which are smooth on scales of 
order 100ft. -1 kpc, their origin and their evolution are largely 
unclear. Fields of that scale must be related to the formation 
and evolution of the host cluster rather than the kinematics and 
internal evolution of the cluster galaxies. It is therefore an in- 
teresting question in what way galaxy clusters deform, com- 
press and amplify the seed fields pervading the material that 
finally ends up in cluster cores, what magnetic field configura- 
tions are prevalent in final, virialised clusters, and whether the 
cosmological background model has any significant effect on 
the evolution of intracluster fields. 

We address these questions with cosmological, magneto- 
hydrodynamic cluster simulations. The numerical code used 
for this study uses smooth-particle hydrodynamics (SPH) for 
the magnetised intracluster gas and purely gravitational inter- 
action between the particles of the dominant dark-matter com- 
ponent. The technical and numerical properties of our "Grape- 
MSPH" code, its performance in various tests, and results of 
a first application to cluster evolution were described in an 
earlier paper (Dolag et al. 1999). Here, we extend this earlier 
study in several ways. First, we perform cluster simulations in 
two different cosmological models, a standard-CDM (SCDM) 
model with ft = 1 and A = 0, and a ACDM model with 
£1 = 0.3 and A = 0.7, to explore the influence of the cosmo- 
logical background model. Second, we focus on the structure 



Send offprint requests to: K. Dolag 



of the intracluster magnetic fields, quantified by radial field- 
strength profiles, correlation and field-reversal lengths, and 
the magnetic-field power spectrum. Third, we investigate how 
magnetic fields grow during cluster evolution, and how they are 
affected by cluster mergers. Finally, we compare the Faraday 
rotation predicted by our cluster models with the largest avail- 
able sample of Faraday-rotation measurements. 

2. Numerical technique and initial conditions 

2. 1. GrapeMSPH 

Our cosm ological MHD code "Grape-MSPH" was described 
in detail in Dolag et al. (1999) Therefore, we only give a brief 
summary here. The code combines the merely gravitational in- 
teraction of a dominant dark-matter component with the hy- 
drodynamics of a magnetised gaseous component. The grav- 
itational interaction of the particles is evaluated on GRAPE 
boards (cf. Sugimoto et al. 1990), while the gas dynamics is 
computed in the SPH approximation (Lucy 1977; Monaghan 
1992). The original "GrapeSPH" code (Steinmetz 1996) was 
augmented by an implementation of the (ideal) magneto- 
hydrodynamic equations describing the evolution of the mag- 
netic fields carried by the gas. Throughout, we assume that 
the electric conductivity of the gas is infinite, which implies 
that the magnetic fields are frozen into the gas flow. The back- 
reaction of the magnetic fields on the gas via the Lorentz force 
is fully included. The numerical viscosity required by SPH to 
properly capture shocks is chosen such that angular-momentum 
transport in presence of shear flows is carefully controlled. The 
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Table 1. Parameters of the SCDM and the ACDM model. 
The SCDM clusters are COBE normalised, which is in con- 
flict with the local cluster abundance. However, this does not 
affect our specific results because we are studying the mag- 
netic field growth in individual clusters. Initial conditions are 
kindly provided by: (1) Bartelmann & Steinmetz 1996, (2) Eke 
et al. 1998. 



SPH kernel width is automatically adapted to the local number 
density of SPH particles, which results in an adaptive spatial 
resolution of the code. 

Extensive tests of the code were performed and described in 
Dolag et al. (1999) The code succeeds in solving the co-planar 
MHD Riemann problem posed by Brio & Wu (1988). Although 
the simulated magnetic field is not strictly divergence-free, 
monitoring shows that V • B always remains negligible com- 
pared to t he magnetic field d ivided by a typical length scale 
of B (see frjolag et al. (1999) for more detail). The code also 
assumes the intracluster medium to be an ideal gas with an 
adiabatic index of 7 = 5/3. Although foreseen and imple- 
mented in the code, we neglect cooling here. One reason is that 
we cannot avoid the cooling catastrophe without a substantial 
source of heat. Most of the gas cools and ends up in an unphys- 
ically condensed central lump. Moreover, the cool gas fraction 
strongly depends on the resolution of the simulation. In addi- 
tion, a substantial amount of numerical heating is expected due 
to our limited mass resolution. 

The surroundings of the clusters are important because their 
gravitational tidal fields affect the overall cluster structure and 
the merger history of the clusters during formation. Therefore, 
the cluster simulation volumes are surrounded by a layer of 
boundary particles whose purpose it is to accurately represent 
the tidal fields of the cluster neighbourhood. 
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2.2. Initial conditions 

We set up initial conditions for two CDM-dominated model 
universes (standard CDM, SCDM; and flat, low-density CDM, 
ACDM), whose parameters are listed in Table 0. For the SCDM 
models, the central region of ~ 70 Mpc (comoving) are filled 
with dark-matter particles with a mass of 3.2 x 10 11 Mq each, 
mixed with an equal number of gas particles whose mass is 
twenty times smaller. The central region is surrounded by col- 
lisionless boundary particles whose mass increases outward to 
mimic the tidal forces of the neighbouring large-scale struc- 
ture. Including the region filled with boundary particles, the 
simulation volume is a sphere with a (comoving) diameter of 
~ 300 Mpc. For the ACDM models, the dark matter particles 
have masses of 8.1 x 10 10 M© each, while the mass of the gas 
particles is ten times smaller. Including the region filled with 
boundary particles, the simulation volume in this case has a 
comoving diameter of ~ 600 Mpc. The simulations labelled 
as "double" are taken from the SCDM models, but have twice 



Table 2. Initial magnetic fields (Col. 2) and the resulting final 
mean magnetic field strengths in the cluster cores (Cols. 3 and 4 
for SCDM and ACDM models, respectively). The final values 
are an average over the gas particles central region (enclosed 
by a radius of 0.1 rw) across the sample of ten clusters simu- 
lated for each cosmology. Note that in the case of doubled mass 
resolution, the amplification of the initial field is significantly 
increased because shear flows are better resolved. This means 
that the initial field required to reproduce the observations with 
our simulations are only upper limits. Apart from that, there is 
no indication that the increased resolution changes the charac- 
teristics of the final magnetic field. 



the number of particles with half the original mass. For each 
cosmology, we used ten different realisations of the density- 
fluctuation field at an initial redshift z; n i = 15 (SCDM, kindly 
provided by M. Steinmetz) and z ini = 20 (ACDM, kindly pro- 
vided by J. F. Navarro). Note that the growth of fluctuations 
is slowed down in a ACDM model compared to the SCDM 
model. The different initial redshifts compensate for that, so 
that the initial fluctuations in linear expectation grow by the 
same amount in both cosmologies. 

Each of these realisations evolves such that it contains clus- 
ters of different final masses and in different dynamical states 
at the final redshift z = 0. We simulate each of these clusters 
with up to five different initial configurations of the magnetic 
field, which are listed in Table. 2. In addition, one set (SCDM) 
is computed with its mass resolution increased by a factor of 
two to test for resolution effects in our results. In total, we sim- 
ulate and study 100 cluster models. 

Lacking any detailed knowledge on the origin of primordial 
magnetic seed fields, we explore two extreme cases of initial 
field configurations. In one case {"homogeneous"), we assume 
that the field is initially constant throughout the cluster volume. 
In the other case {"chaotic"), we let the initial field orientation 
vary randomly from place to place, subject only to the condi- 
tion that V ■ B = 0. The initial field strengths in both cases are 
determined by setting the mean field energy densities. To study 
the effects of changing the mean initial energy density of the 
magnetic seed field, we ran simulations using three choices for 
this quantity, labelled "low", "medium" and "high". 

To gain some insight into performance of the numerics, 
we also ran a set of simulations with doubled mass resolution, 
labelled "double". The magnetic field configurations of these 
runs are the same as in the "homogeneous" case. Of course, 
this is not meant to replace a full resolution study. Here, we 
are limited by the capacity of our GrapeMSPH code, because 
the Grape design restricts the number of gas particles. We also 
ran a set of control simulations for each cosmological model 
without magnetic field, labelled "no". In total, we have a set 
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Table 3. Parameters (amplitude A and slope a) of power 
laws best fitting the relation between final mean magnetic 
field strength and mean temperature in the simulated clusters. 
Numbers for different cosmologies, different magnetic field se- 
tups and different radii for the averaging volume are deter- 
mined. The different slopes fitted to the different models are 
all consistent within the model uncertainties. 



of ten realisations for each cosmology, simulated with six dif- 
ferent magnetic field setups in the SCDM cosmology, and with 
four different magnetic field setups in the ACDM cosmology, 
leading to the 100 simulated clusters mentioned before. 

Table ^ summarises the initial field set-ups and the mean 
field strengths in the cluster cores. The initial field strengths are 
of order 10 -9 G, the final field strengths of micro-Gauss order 
are averaged over the gas particles within spheres of 0.1 r v i r 
and across the sample of ten clusters. For some aspects of the 
analysis of the simulations to be described later on, we also use 
ten less massive objects identified close to the main clusters. As 
they are relatively more poorly resolved in the simulations, the 
results from this smaller objects have to be treated with caution. 

Furthermore, we point out that as we are including the 
back-reaction of the magnetic field on the gas, different ini- 
tial field strengths and configurations are changes to the initial 
conditions from which the clusters develop, hence their evolu- 
tion proceeds in a somewhat different way. This implies that 
even a low magnetic field can alter the dynamical time scale 
of a simulated cluster to some degree. Therefore, when com- 
paring simulated clusters with different initial magnetic fields 
at equal redshifts, one has to keep in mind that they can be 
in slightly different evolutionary stages according to their indi- 
vidual dynamical time scale. Of course this has an even more 
pronounced effect on the models with doubled mass resolu- 
tion. In this case, it is not the magnetic field but the addition 
of new particles at intermediate positions which changes the 
initial conditions and leads to slightly different clusters which 
evolve in a slightly different fashion. Figure |l]provides a visual 
impression of the differences between six different simulations 
of one cluster in a pre-merger phase, starting from the same 
SCDM cluster initial conditions. The simulations use four dif- 
ferent configurations of the initial magnetic field, one run with- 
out magnetic field, and one with doubled mass resolution. 
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Table 4. Cols. 2 and 3 show typical field-reversal scales in kpc 
(physical units) for the magnetic fields and their rms deviation 
for the clusters in the SCDM model. The values are calculated 
using all particles in cluster-centric spheres with one tenth of 
the viral radius, and averaged across the sample of ten clusters 
per cosmological model. 



3. Structure of the final magnetic field 

As reported earlier (Dolag et al. 1999), the process of cluster 
formation and the surrounding large scale structure entirely de- 
termine the magnetic field properties in the simulated clusters. 
As shown before, all information contained in the structure of 
the seed field is completely wiped out in the process of the 
collapse. This means that the final magnetic field properties 
obtained reflect our understanding of the formation of galaxy 
clusters, with the underlying general assumption that the mag- 
netic field in galaxy clusters results from the amplification of 
weak seed fields. 

3. 1. Amplification by shear flows 

Apart from the compression of the magnetic field frozen into 
the gas, the magnetic field is also substantially amplified by 
shear flows as the gas is accreted by the cluster. A necessary 
condition for this amplification mechanism is that the velocity 
difference across the boundary layer exceeds the local Alfven 
speed, which is proportional to the field strength, so that the 
condition is more easily satisfied for weaker fields. This means 
that the relative field amplification becomes smaller as the 
fields grow. As the shear flows are largely driven by gradients in 
the dark matter potential, this mechanism is quite independent 
of pressure support by the magnetic field and therefore reduces 
the relative amplification much earlier than local equipartition 
is reached. 

3.2. Final field strength 

One quantity of primary theoretical interest is the final mag- 
netic field strength in galaxy clusters predicted by our mod- 
els. Of course, it depends on the strengths of the seed fields 
assumed in the simulations. The final magnetic field strengths 
also depend on the final cluster mass and on the dynamical state 
of the galaxy cluster because mergers are known to amplify the 
magnetic field (Roettiger et al. 1999, Dolag et al. 1999). 

The simulations also predict that the magnetic field de- 
creases with increasing cluster-centric distance. Therefore, we 
compute the mean magnetic field strengths within 10 per cent 
and 100 per cent of the virial radius of the final cluster. Figure^ 
shows how the mean magnetic field strength scales with tem- 




Fig. 1. The figure shows versions of the same SCDM galaxy cluster for six different initial magnetic field models as indicated 
in the plot, including one version without magnetic field for reference. Plotted are the projected rotation measure maps, an 
iso-density surface and a cut through the magnetic field. 
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Fig. 2. The mean magnetic field is shown as function of the mean mass weighted temperature inside the clusters. Here, the ten 
less massive objects identified in the simulations are included in the analysis. The symbols belong to the "medium" model, the 
lines represent the power laws fitted to all models for the different magnetic seed fields. Both quantities are calculated within 
10 per cent (thin symbols and lines) and 100 per cent (thick symbols and lines) of the virial radius. The left side represents the 
SCDM models, the right side the LCDM models, respectively. 



perature for three sets of initial field strengths in our SCDM 
simulations. In order to increase the temperature range of our 
cluster set, we also include ten less massive objects found near 
the most massive clusters in our simulations. We thus ended up 
with a total of twenty different objects for every initial mag- 
netic field configuration in both cosmologies. 

We fitted this relation in all sets of simulations by a power 
law, (B) = AT a . The values of A and a we obtained for all 
our simulations in both cosmologies are summarised in Table |. 
In Fig. |5[ the different symbols and line types distinguish be- 
tween the different initial field models. The thin and thick sym- 
bols and lines show results calculated within 10 per cent and 
100 per cent of the virial radius, respectively. 

As Table || shows, this relation tends to flatten when the 
magnetic field strength is averaged within a larger volume, and 
also when the initial seed magnetic field strength is increased. 
These trends are not significant given the uncertainties of the 
models, but they may have the same reason as the flattening 
of the magnetic field profiles towards the cluster cores, which 
will be discussed below. It is simply the fact that smaller mag- 
netic fields can more easily be amplified by shear flows than 
already grown magnetic fields, as discussed before. These rela- 
tions again show that there is no difference in the final states 
of the simulation sets with "chaotic and "homogeneous ini- 
tial field configurations, which confirms that the information 
on the initial field configuration gets entirely wiped out during 
the process of cluster collapse. There is also no visible differ- 
ence between the two cosmologies. Doubling the mass resolu- 



tion does not significantly change the final field structure either, 
although the field amplification is higher. 

Another interesting result of the simulations is the relation 
between the magnetic field strength and the gas density in the 
clusters. Figure ^ shows the mean magnetic field strength in 
spherical shells centred on the cluster centre as function of 
cluster-centric distance in units of the virial radius. Plotted are 
the profiles at redshift zero, but the shape does not change sig- 
nificantly between redshift one and zero. The gas density pro- 
files, averaged across the sample of ten clusters and raised to 
the power of 2/3, are over-plotted in arbitrary units to match the 
magnetic field strength amplitude at the virial radius (dashed 
lines). 

The average magnetic field strength profiles in the outer 
regions of the clusters (beyond ~ 0.3r V i r ) follow the density 
profile very well, while they flatten off towards the cluster core, 
where the magnetic field profiles are shallower than the density 
profiles. Plotting the density profile raised to the power of 2/3 
demonstrates that in the outer parts of the cluster, the average 
magnetic field strength is more strongly amplified as expected 
from compression alone (oc p 2 / 3 ), thus shear flows have to play 
an important role there. In the centre, the field amplification is 
reduced. This is due to the fact that the higher magnetic field 
strength in the cluster centres slows down the further field am- 
plification compared to the outer cluster regions. There are two 
reasons for that. First, the motions of the gas are less coordi- 
nated in the cluster cores compared to the large-scale motion 
of the infalling gas in the outer parts, so that amplification by 
ordered shear flows is suppressed. Moreover, increasingly im- 
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Fig. 3. The two panels show the radial profiles of the magnetic field strengths (thick lines) in the simulated, final clusters. The 
profiles are obtained by averaging across the sample of ten clusters. The averaged gas density profiles, raised to the power of 
2/3, are over-plotted for comparison as the thin lines, arbitrarily normalised to the amplitude of the magnetic-field profiles at the 
virial radius. The left panel shows the SCDM simulations for the various initial field configurations, as labelled in the plot. The 
right panel shows the results for the ACDM models. 



portant back-reaction of the field on the gas flow counteracts 
further amplification. Second, the magnetic field is also sub- 
stantially amplified by shear flows as the gas is accreted by the 
cluster. A necessary condition for such an amplification is that 
the relative velocity on both sides of the boundary layer ex- 
ceeds the local Alfven speed which is proportional to the field 
strength, so that the condition is more easily satisfied for the 
weaker fields embedded in the infalling material in the outer 
regions of the clusters than for the stronger fields near cluster 
centres. 

This is also responsible for flattening the relation between 
temperature and magnetic field strengths for increasing ini- 
tial magnetic field strength. Table || lists the final magnetic 
field strengths in the cores of the simulated cluster sample for 
the different initial magnetic field configurations. Again, the 
chaotic and the homogeneous initial field configurations lead 
to the same radial magnetic field strength profiles. As observed 
before in Dolag et al. (1999), the final magnetic field struc- 
ture is entirely driven by the collapse of the galaxy cluster, and 
any information about the initial field structure is completely 
erased. This justifies the use of homogeneous initial conditions 
for the magnetic field setup, as they are numerically easier to 
handle. The profiles do not change significantly when the par- 
ticle number is doubled. 

These results on the radial magnetic field strength profiles 
and the dependence of mean field strengths on the cluster tem- 
perature may help to reconcile the observations in different 
galaxy clusters, which seem to give somewhat different val- 
ues for the magnetic field for different clusters observed. In 



addition, the difference between magnetic field inferred from 
Faraday rotation measurements on the one hand (testing the 
line-of-sight average) and the combination of radio emission 
with hard X-ray observations on the other hand (measuring vol- 
ume averaged magnetic fields) could be bridged using our sim- 
ulations, although the relativistic electron density needs to be 
assumed. Comparing the mean magnetic field within the virial 
radius and one tenth of the virial radius, the data points in Fig. |2| 
show that these values differ already by factors of five, or more. 

3.3. Field reversals 

The formation of large scale structure in the surroundings of 
a cluster also imprint some structure on the magnetic field of 
galaxy clusters. Large scale flows can order magnetic fields 
parallel with the flow. Also the field compression in shock 
fronts orders the magnetic fields to some degree. Since the 
shocks compress the medium in only one direction, the mag- 
netic fields in the shocked gas tend to lie within the surface of 
the shock shell. On the other hand, the more chaotic motion of 
the gas inside cluster cores tends to destroy the order imposed 
on the magnetic field during cluster collapse. We thus expect 
a complicated behaviour of the field structure and a range of 
autocorrelation length scales. It is of interest for the interpre- 
tation of rotation measures on what scales field orientations 
reverse inside cluster cores. Our simulations predict how mag- 
netic fields can be ordered in the process of large scale structure 
formation and gravitational collapse. In addition, there might 
be smaller scale contributions from the outflows and the co- 
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Fig. 4. The figure shows the autocorrelation function of the 
magnetic field, (B(x) ■ B(x + d)) as a function of \d\ for 
one SCDM cluster, simulated with five different initial condi- 
tions for the magnetic field, as indicated in the plot. The num- 
bers attached to the minima and maxima indicate the number of 
particle pairs used to compute the mean value. Only particles 
within ten per cent of the virial radius are taken into account. 
All curves are normalised to fall within ±1 within the region 
shown. Table |] lists the mean inferred autocorrelation length 
scale across our cluster sample and its rms deviation. 

coons of the radio galaxies, from which Faraday rotations are 
measured. 

In order to infer this typical length scale in the simulated 
clusters, we cut out the central sphere with one tenth of the 
virial radius of each cluster and calculate the autocorrelation 
function 

((d) = (B(x)-B(x + d)) (1) 

by identifying all possible pairs of particles, sorting them by 
their separation, and binning them. The mean is then calculated 
over the number of particles in each distance bin. This function 
should be positive for scales smaller than that of typical field re- 
versals, then go through a minimum at this typical length scale, 
and then damp due to the mixing of different length scales. This 
mixing, the presence of shocks, and the layers between shear 
flows, already affect this measurement at the length scale of in- 
terest, which render the determination (and even the meaning) 
of the length scale dubious. On larger scales, it gets confused 
by the mixture of different length scales present in the cluster. 
Therefore, the minimum correlation length scale is difficult to 
measure. 

Figure ^ shows the autocorrelation function for one simu- 
lated SCDM cluster. At small scales, this function shows the 
positive value expected, but it is already very noisy. Note that 



Fig. 5. This figure shows the magnetic field component in the x 
(solid line) and y directions (dashed line) along a line-of-sight 
through one of the simulated clusters. This line was chosen par- 
allel to the x axis and passing through the cluster centre. The 
envelope of the curves follows the general radial decrease in 
the magnetic field strength. It is clearly seen that field reversals 
are relatively more frequent in the cluster centre, and relatively 
less frequent away from the cluster core. 



the innermost distance bin is narrower than the average inter- 
particle distance, and are therefore likely to be in shocks. 

Although the magnetic field in the shocks is preferentially 
aligned tangentially to its surface, the orientation within this 
surface can be very tangled. Accordingly, the innermost bin is 
very noisy and can even yield negative values, as seen in the 
"low" model. At somewhat larger distances, the autocorrelation 
function becomes positive, as expected. Then, at almost twice 
this distance, £(d) becomes negative, also as expected. Clearly, 
the field autocorrelation function is determined by the structure 
of the cluster and shows a very similar behaviour for differ- 
ent initial magnetic field configurations in that it consistently 
reaches a pronounced minimum between 60 and 80 kpc. The 
numbers attached to the individual lines denote the number of 
particles used to calculate the mean within this distance bin. At 
the first minimum, several hundred particles are used to com- 
pute the SPH average. The autocorrelation function reaches 
zero only beyond the plotted range. 

Since the autocorrelation signal is the result of mixing 
structures with different length scales (due to substructure), re- 
gions simulated at varying spatial resolutions (because of the 
adaptive SPH kernel), and volumes with varying magnetic field 
strength (due to the field structure in the cluster), it is not at all 
straightforward to read off a characteristic autocorrelation scale 
from such plots for all clusters. In order to obtain a more robust 
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estimator, we tried to read off the first two extrema and the zero, 
and constructed the appropriately weighted mean 

I = (2 /drop + /zeroj + 2/ 3 / m iu + 1/2 Z ZC ro 2 )/4 (2) 

of the positions of the extrema to get the length scale / rather 
than trying to identify the location of the first minimum only. 
Interpreting the signal like a wave and taking the distance be- 
tween the zero crossing as the length-scale of the magnetic 
field, Zdrop in this equation is the length scale where the sig- 
nal starts to decrease, which is an indicator of half the typical 
length of the magnetic field / zorol , where the signal is expected 
to drop to zero. The next minimum / m ; n should reflect 2/3 of 
the magnetic field reversal scale. We may also get a second 
minimum, Z zcr02 , reflecting twice the magnetic field reversal 
scale, before the signal fades out due to the reasons mentioned 
above. 

Table [f] summarises the field reversal scales found in the 
simulations. The table lists the averages across the sample of 
ten clusters for each set of simulations and their rms deviations. 
Clearly, the reversal scale does not depend on the structure and 
the strength of the initial magnetic field. Doubling the particle 
resolution does not yield significantly different answers either. 

3.4. Field structure 

The length scale on which the magnetic field typically reverses 
its direction is important for the interpretation of Faraday ro- 
tation measurements. Since the Faraday rotation measure is an 
integral of the magnetic field component parallel to the line- 
of-sight, weighted by the electron density, it is also important 
to know how the field reversal scale changes across the cluster. 
Figure || shows one line of sight passing through the centre of 
one of the simulated clusters. The solid and dashed lines show 
the x component (along the line-of-sight) and y component 
(perpendicular to the line-of-sight) of the magnetic field along 
the line-of-sight, respectively. The line-of-sight was chosen to 
go through the centre of mass of the cluster and follows the x 
axis. Notice that the plot is composed of four panels, where the 
regions near zero are cut out to allow a logarithmic scale. 

It can be read off the figure how frequently the magnetic 
field changes it direction, and how the reversal length scale 
changes with the distance from the cluster centre. The envelope 
of the curves in this plot nicely follows the radial decline of the 
magnetic field strength in our simulations. It is also clearly seen 
that the magnetic field reverses relatively less often in the outer 
parts of the cluster. While we count 23 extrema for the x com- 
ponent of the magnetic field within ±4 Mpc (in physical units), 
there are already 15 extrema within ±1 Mpc (in physical units). 

This effective increase in the reversal length scale has to be 
taken into account equally well as the decrease of the magnetic 
field strength with radius when inferring magnetic fields from 
Faraday rotation measurements. These effects tend to reduce 
the magnetic field strengths inferred from current observations. 

In order to improve our understanding of the statistical 
properties of the field structure, we calculated the power spec- 
tra of the magnetic field in the simulated galaxy clusters. To 
this end, we calculated the auto-correlation function of the field 
strength within a cube of 10 Mpc comoving side length, centred 
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Fig. 6. This figure shows the magnetic-field power spectra cal- 
culated from the simulation of one galaxy cluster for different 
initial magnetic field configurations, distinguished by different 
line types as indicated in the plot. For each simulation, the best- 
fitting power law is overlaid. 



model 


a 


b 


(b) 


RMS 6 


low 


9.9 


-2.7 






chaotic 


10.6 


-2.7 






double 


8.5 


-3.0 






medium 


6.8 


-2.7 


-2.7 


0.2 


high 


3.5 


-2.4 







Table 5. Best-fitting parameters for power laws, a k b , to the 
magnetic-field power spectra. The table summarises the values 
for different initial magnetic field configurations in the SCDM 
models. 



on the cluster centre, using an equidistant grid of 80 3 cells on 
which all quantities are calculated by means of the SPH for- 
malism. 

In this type of simulation, the dynamical range of spatial 
scales is not high enough to fully resolve the correlation func- 
tion over a satisfactory range of length scales. Nonetheless, we 
could compute the correlation function and the power spectra 
over one and a half orders of magnitude in length scale, within 
which it nicely follows a power law. The hump towards small 
k values (i.e. large radii) can qualitatively be understood. First, 
the correlation function rises towards the typical distance be- 
tween major objects within the simulation, e.g. the cluster and 
the next major object falling onto it. In addition, the correlation 
function in a cubical simulation volume is over-predicted for 
scales approaching half the cube side length (i.e. the Nyquist 
length) because one particle of each pair is most probably po- 
sitioned near the cluster centre. On even larger scales, one par- 
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ticle of each pair most likely falls outside the cluster core, 
and therefore the correlation function drops rapidly. On small 
scales, the resolution of the cube used, and the resolution of 
the simulation itself, restrict the calculation of the correlation 
function and thus also the computation of the power spectra. 

Nonetheless, it is possible to calculate the power spectra 
over a reasonably wide range in wave number, and compare 
the results for different magnetic field configurations. Figure || 
shows the power spectrum calculated from one of the simulated 
clusters for different initial magnetic field configurations. The 
power spectra look quite the same for all models. 

The exponents of the best-fitting power laws for one set of 
simulations range from —2.3 to —3.1, with a mean of —2.7. 
This value is surprisingly low. For Kolmogorov turbulence, 
—5/3 would be expected. A fully established MHD turbulence 
would produce —3/2, and even in the extreme case of a young, 
incomp letely established MHD tu rbulence, we would expect 
—2 (cf. Eilek & Henriksen (1984)| and references therein). 

There are two cases discussed in the literature which result 
in a smaller power- law slope. First, f or a 2-dimensional Navier- 
Stokes turbulence ( Biskamp (1993)| ), an exponent of —3 would 
be expected. As shear flows play a major role in amplifying the 
magnetic field, it is very likely that we see our cluster simula- 
tions in a state, where the turbulence induced by shear flows 
evolves into a Kolmogorov turbulence. The second case, which 
leads to a power-law index of —4 (Eilek & Henriksen (1984)), 
is turbulence in a very viscous medium. SPH is known to be 
viscous, so this may be a possibility to explain the small expo- 
nent found, but as the artificial viscosity in SPH is constructed 
to vanish in shear flows, this is unlikely to be the complete ex- 
planation. To clarify this issue, detailed tests remain to be done. 
Table || summarises the power-law fits to the power spectra in 
the different models. In particular, the spectra for the chaotic 
and the equivalent homogeneous initial field configurations are 
indistinguishable. Also, the power spectra for the run with dou- 
bled mass resolution do not differ, which argues against a dom- 
inant impact of the SPH viscosity on the behaviour of the mag- 
netic field in the simulations. There is also no visible change in 
the power laws with changes in the magnetic field strength. 



4. Magnetic field evolution 

This section describes how the magnetic fields evolve during 
the formation of galaxy clusters in our simulations. The field 
evolution is driven by structure formation; in particular, major 
mergers have a strong influence on the how the field evolves. 

4. 1. Magnetic field evolution 

We follow two different approaches for describing the evolu- 
tion of the magnetic field in our simulations. First, we look 
at each individual given cluster and measure how the mag- 
netic field strength in the cluster (e.g. within the core) evolves 
with time. On the other hand, we can ask how the magnetic 
field evolves which is carried by the gas that finally ends up 
in the cluster core. The first approach could be considered as 
Eulerian, the second as Lagrangian. 
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Fig. 7. The mean magnetic field in one cluster simulated in the 
SCDM cosmology (h = 0.5) of those gas particles which at the 
end of the simulation populate spheres of ten and fifty per cent 
of the virial radius (dotted and dashed lines) around the cluster 
centre is shown as a function of redshift. The solid line is the 
mean magnetic field within a sphere of comoving radius 1 Mpc 
which remains centred on the most massive cluster component. 

Figure ^ shows the results of both approaches for one clus- 
ter simulated in the SCDM cosmology (h = 0.5)). The solid 
line shows the mean magnetic field within a sphere of 1 Mpc 
radius around the cluster centre. More precisely, it shows the 
mean magnetic field within a sphere of 1 Mpc comoving ra- 
dius, centred on the most massive halo contributing to the final 
cluster. The dotted and dashed lines show as functions of time 
the mean magnetic fields carried by all particles which finally 
fall into spheres around the cluster centre with radii of ten and 
fifty per cent of the virial radius. It is clearly seen that the mag- 
netic field in the core of the evolving galaxy cluster is stronger 
than the magnetic field in the cluster surroundings, which are 
finally accreted by the galaxy cluster. All three curves show 
pronounced maxima which coincide with major mergers at red- 
shifts 0.5 and 0.2, and likewise at the second core passage of 
a major merger at redshift 0.7. Also, the combined effect of a 
second core passage and a merger at redshift ~ 0.2 can be seen 
here, and more explicitly in Fig. |8| The magnetic field in the 
cluster cores and in related shocks during mergers is increased 
substantially. Even the mean field within a large volume can be 
increased by a factor of a few, as can be read off Fig. ^| 

4.2. Influence of mergers on observables 

We now investigate how the projected patterns of Faraday ro- 
tation measures evolves. In doing so, we calculated synthetic 
rotation measure maps for 60 epochs between redshifts 0.8 and 
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SCDM 


5 


20 


50 


low B 


-0.23 


-0.33 


-0.31 


medium B 


-0.09 


-0.15 


-0.15 


highB 


-0.04 


-0.07 


-0.08 


ACDM 


5 


20 


50 


low B 


-0.29 


-0.34 


-0.42 


medium B 


-0.16 


-0.18 


-0.24 


highB 


-0.09 


-0.11 


-0.15 



Table 6. This table summarises the slopes of the linear 
fits shown in Fig. ^ for different cosmologies and different 
Faraday-rotation thresholds. 



(today). For each epoch, we have one map for each of the 
three different projection directions (chosen to be the coordi- 
nate axes). In each map, we calculate the total area of those 
regions in which the rotation measure exceeds a given thresh- 
old (in other words, we compute the measure of the excursion 
set of the Faraday rotation). 

Figure ^| shows this area as a function of time (or redshift), 
normalised to the area which emits 95% of the X-ray lumi- 
nosity. The different line types distinguish between the three 
different projection directions, the three bundles of lines are 
for rotation measure thresholds of 5, 20, and 50 rad/m 2 . The 
general tendency for this quantity to grow with time is evident, 
as is expected from the growing magnetic field. The structure 
seen in the time evolution is clearly related to merger events. 
All peaks on these curves correspond to a core passage dur- 
ing a merger event. Likewise, the second core passage after the 
merger causes the area to peak. The signal does not strongly 
depend on the projection direction, and the sequence of peaks 
is present for all rotation-measure thresholds chosen. This pro- 
vides clear evidence that merger events strongly contribute to 
the amplification of the magnetic field, an effect which is con- 
firmed by other simulations (Roettiger et al. 1999). Since this 
effect is visible in all projections, this also applies to mergers 
along the line-of-sight, whose X-ray contours look quite regu- 
lar. This implies that determining the area covered by Faraday 
rotation measures above a fixed threshold may help constrain- 
ing cluster mergers in such configurations. 

4.3. Evolution in different cosmologies 

Figure ^ shows the same area, but averaged over all ten real- 
izations. The sequence of peaks is smoothed out, leaving the 
general trend. The different initial field configurations are indi- 
cated by different types of symbol, as indicated in the plot. We 
overlay fits showing that the growth is very well fit by linear 
functions. Table ^summarises the gradient inferred from these 
fits for different cosmologies, different initial field configura- 
tions and different thresholds. 

5. Comparison with measurements 

Figure [l(] compares the measurements of Clarke et al. (2001) 
with the synthetic Faraday rotation measures taken from 
our simulations. These new measurements are in excellent 
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Fig. 8. This figure shows how the area, covered by a Faraday 
rotation measure exceeding thresholds of 5 (upper bundle), 20 
(middle bundle) and 50 rad/m 2 (lower bundle) evolves with 
time. The area is normalised by the area emitting 95% of the 
cluster X-ray luminosity. The different line types distinguish 
the three different projection directions. 

agreement with previous measurements collected by Kim et 
al. (1990), and are also compatible with Faraday-rotation mea- 
surements in individual clusters like Abell 119 by Feretti et 
al. (1999). 

To calculate the synthetic signal, we computed three 
rotation-measure maps for the three spatial projects for each 
cluster on a grid of 60 2 cells, covering 5 x 5Mpc (physical 
units). For each cosmology and all initial magnetic field con- 
figurations, the Faraday-rotation measures are sorted in radial 
bins and combined over all ten realizations. The simulated val- 
ues compare very well with the observations. We calculated the 
signal from the observations between 3 Mpc and 10 Mpc as a 
reference (11 rad/m 2 ) and added this as a background to the 
synthetic data. The synthetic curves very well follow the shape 
of the measurement distributions. It is also easily seen that the 
initial magnetic field strength required to reproduce the obser- 
vations falls somewhere between the "medium" and the "high" 
models. 

Again, it is obvious that the difference for the two different 
initial field configurations, namely homogeneous and chaotic, 
is negligible. Also, the two different resolutions lead to sta- 
tistically the same signal, but as the mass resolution is only 
changed by a factor of two, this does not resemble a conver- 
gence study. 

6. Summary and conclusions 

We performed cosmological, magneto-hydrodynamic simula- 
tions of galaxy clusters to study the evolution and final struc- 
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Fig. 9. This figure shows the same quantity as shown in Fig. §, but only for the one threshold of 20 rad/m 2 , and averaged over the 
three projection directions and over the ten clusters per sample. The different line types represent the different sets of models for 
the initial magnetic field. The left panel shows results for the SCDM cosmology, the right panel for the ACDM cosmology. 



ture of magnetic fields frozen into the intracluster medium. In 
total, we investigated 100 cluster models of different masses 
in two different cosmologies (standard and open CDM), vary- 
ing the structure and the strength of the initial magnetic field 
configuration. Our results can be summarised as follows: 

- Starting with magnetic fields of order 10~ 9 G at an initial 
redshift of z ini = 15 (20 for ACDM), final fields reach 
micro-Gauss strength in cluster cores at redshift z — 0. 
This amplification factor of ~ 10 3 was found and discussed 
earlier (Dolag et al. 1999). It exceeds the amplification fac- 
tor expected from spherical collapse by an order of mag- 
nitude, which can be traced back to field amplification in 
shear flows. 

- Radial profiles of the magnetic field strength closely follow 
the cluster density profiles at radii of ~ 200h~ 1 kpc and 
beyond. At smaller radii, the magnetic field strength pro- 
files flatten. This reflects the behaviour of the gas distribu- 
tion, which also develops a flat core despite the dark-matter 
density cusp because of the finite gas pressure. In a purely 
spherical collapse, the magnetic field strength is expected 
to scale with density raised to the 2/3 power, flattening 
the profile. Additionally, the amplification by shear flows 
works more efficiently when the field strengths are small 
and therefore the fields get more relative amplification in 
the outer parts, leading to further flattening of the magnetic 
field profile compared to the gas profile. This happens in- 
dependently of the strength and the structure of the initial 
magnetic seed field, and for both cosmologies investigated. 

- The correlation lengths of the magnetic fields in cluster 
cores are of order ~ 50ft- -1 kpc. Similarly, the length scale 
of field reversals along lines-of-sight through cluster cores 



is of order ~ 100ft -1 kpc. These length scales show no 
significant dependence on the initial field strength or con- 
figuration, or on cosmology. 

- The power spectrum of the magnetic field can closely be 
approximated by a power law, Pb (k) oc k~ n , with a power- 
law exponent of n ~ —2.7. This is steeper than expected 
for Kolmogorov or MHD turbulence, but close to the ex- 
pectation for two-dimensional Navier-Stokes turbulence. A 
possible explanation is the occurrence of essentially two- 
dimensional shear flows and their importance for the mag- 
netic field evolution in clusters. Although we do not expect 
this to be a strong effect, the high numerical viscosity in 
SPH codes also tends to steepen the power spectrum. 

- Predominantly driven by merger events, the mean magnetic 
field in cluster cores grows roughly exponentially with red- 
shift. Between redshifts z = 0.8 and z = 0, the growth is 
well approximated by 

\B\(z) w 10 -2 - 5z mG . (3) 

Superposed on this general trend is the effect of mergers, 
during which the magnetic field strength goes through pro- 
nounced maxima. 

- Field amplification by mergers shows up prominently in ob- 
servable quantities like the Faraday rotation. During merg- 
ers, the projected cluster area in which substantial Faraday 
rotation occurs rises sharply, and it drops during the subse- 
quent relaxation phases. 

- Simulated Faraday rotation measures show that all our clus- 
ter models are in excellent agreement with the largest avail- 
able sample of observed Faraday-rotation measurements, 
although the samples are still too small to provide narrow 
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constraints on, e.g., the initial strength of the magnetic seed 
fields. 

Our results therefore show that the strength of intraclus- 
ter magnetic fields grows rapidly with redshift, and it closely 
follows the cluster density profiles outside a core region of 
~ 200/i _1 kpc radius. The structure of the field can be char- 
acterised by an auto-correlation or reversal scale of order 
50 — lOO/i -1 kpc, well above the spatial resolution of our nu- 
merical simulations. The power spectrum of the magnetic field 
falls steeply, but the distribution of energy density with respect 
to k scales as k 3 P(k) and thus remains monotonically increas- 
ing. This indicates that most of the magnetic energy density is 
located on small scales. 

The independence of these results on cosmology and ini- 
tial field structure demonstrates that the final intracluster fields 
are entirely dominated by the cluster collapse, thus confirming 
and extending the previous findings (Dolag et al. 1999). The 
importance of mergers for observable signatures of intracluster 
magnetic fields, and their steep evolution with redshift, suggest 
important observational tests for the scenario underlying the 
formation and evolution of our simulated galaxy clusters. 
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Fig. 10. This figure shows a comparison of observed Faraday rotation measures with our simulations. Plotted is the absolute 
value of the observed rotation measures as a function of the distance to the next Abell cluster (crosses), taken from Clarke 
et al. (2001). The solid line marks the sample median of the observations. The superposed lines are calculated from synthetic 
Faraday rotation measurements taken from our simulations in an attempt to mimic the observations. The different line types 
distinguish between the different initial magnetic field configurations used in the simulations. The left and right panels show 
SCDM and ACDM results, respectively. The Medians were calculated in bins containing 15 data points each. The shape of 
the synthetic measurements in both cosmologies follow very well the observations. To perfectly match the observations in this 
comparison, a magnetic field somewhere in between the "medium" and the "high" models is needed. 



